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ABSTRACT 

We present new moderate-resolution, far-ultraviolet spectra from the Hubble Space Tele- 
scope/Cosmic Origins Spectrograph (HST/COS) of the BLLac object lES 1553-1-113 covering the 
wavelength range 1135 A < A < 1795 A. The data show a smooth continuum with a wealth of 
narrow (6 < 100 km s~^) absorption features arising in the interstellar medium (ISM) and intergalac- 
tic medium (IGM). These features include 41 Lya absorbers at < ^abs < 0.43, fourteen of which 
are detected in multiple Lyman lines and six of which show absorption in one or more metal lines. 
We analyze a metal-rich triplet {lS.cz ~ 1000 km s^^) of Lya absorbers at Zabs ~ 0.188 in which 
O VI, N V, and C III absorption is detected. Silicon ions (Si III, Si IV) are not detected to fairly 
strong upper limits, and we use the measured Si III/C III upper limit to derive an abundance limit 
(C/Si) > 4 (C/Si)0 for the strongest component of the absorber complex. Galaxy redshift surveys 
show a number of massive galaxies at approximately the same redshift as this absorption complex, 
suggesting that it arises in a large-scale galaxy filament. As one of the brightest extragalactic X-ray 
and 7-ray sources, lES 1553-1-113 is of great interest to the high-energy astrophysics community. With 
no intrinsic emission or absorption features, lES 1553-1-113 has no direct redshift determination. We 
use intervening Lya absorbers to place a direct limit on the redshift: Zom > 0.395 based on a confirmed 
Lya+0 VI absorber and Zcm > 0.433 based on a single-line detection of Lya. The current COS data 
are only sensitive to Lya absorbers at z < 0.47, but we present statistical arguments that z^m ^ 0.58 
(at a Icr confidence limit) based on the non-detection of any Ly/3 absorbers at z > 0.4. 

Subject headings: BL Laccrtae objects: individual: lES 1553-1-113, galaxies: active, intergalactic 
medium, quasars: absorption lines, ultraviolet: general 
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1. INTRODUCTION 

The current interpretation of BLLac objects (Ghis- 
ellini, Maraschi, & Treves 1985) is that they are ac- 
tive galactic nuclei (AGN) with a strongly relativistic jet 
pointed toward our line of sight. As such, any line emis- 
sion or accretion disk features seen in most other types 
of AGN could be masked by the bright jet if present 
in BL Lac objects. Their spectra usually show a fea- 
tureless power-law continuum extending from radio to 
X-ray wavelengths. This spectral characteristic makes 
BL Lac objects ideal for observing intervening absorp- 
tion features arising in the interstellar medium (ISM) 
and intergalactic medium (IGM). Since their continuum 
is easily defined, they make excellent targets for studying 
weak metal-line systems and low-contrast, highly ther- 
mally broadened H I absorbers (e.g., Richter et al. 2004; 
Lehner et al. 2007; Danforth, Stocke & Shull 2010). 

The BLLac object lES 1553-1-113 shows the charac- 
teristic featureless power-law spectrum and is one of the 
brightest known sources of extragalactic high-energy ra- 
diation from X-rays up to VHE (TeV) photons (Costa- 
mante & Ghisellini 2002). However, the featureless 
spectrum makes it difficult to determine the redshift of 
the object and hence its luminosity. Indirect methods 
have given a wide range of limits for the redshift of 
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lES 1553-1-113; the nondetection of a host galaxy gave 
limits from Zgm > 0.09 to z^m > 0.78 (Hutchings & Neff 
1992; Scarpa et al. 2000; Urry et al. 2000; Carangelo et al. 
2003; Sbarufatti et al. 2006; Treves, Falomo, & Uslenghi 
2007). The shape of the 7-ray spectrum observed by 
the Fermi Observatory and ground-based VHE detectors 
(HESS, MAGIC) constrains the redshift to values from 
•Zom < 0.4 to Zom < 0.8 (Aharouiau et al. 2006; Albert 
et al. 2007; Mazin & Goebel 2007; Abdo et al. 2010) 
based on assumptions about the intrinsic spectral energy 
distribution (SED) and pair-production interactions with 
the cosmic infrared background. The only direct redshift 
determination (zem = 0.36; Miller & Green 1983) was 
based on a spurious feature in low-resolution UV spectra 
from the International Ultraviolet Explorer (lUE). The 
detection was later retracted (Falomo & Treves 1990), 
but the erroneous redshift value lives on. 

lES 1553-1-113 is of interest as a bright background 
continuum source for detecting intergalactic absorption 
along the sight line. Bright X-ray sources are especially 
valuable for potentially detecting the long-predicted 
O VII and O VIII tracers (Bregman 2007) of intergalactic 
gas at T = 10^ — 10^ K. Even for a bright X-ray source, 
the required integration times would be very long. How- 
ever, a sufficiently long IGM pathlength provided by a 
bright high-z target would make the required observing 
time investment more attractive. 

In this paper, we present the first medium-resolution 
far-UV spectroscopic observations of lES1553-fll3 in- 
cluding Hubble Space Telescope/ Cosmic Origins Spectro- 
graph (HST/COS Green et al. 2010; Osterman et al. 
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Fig. 1.— Overview of the COS spectrum of lES 1553+113. The full COS/G130M+G160M datasct is shown, smoothed by 35 pixels (~ 5 
resolution elements). Error is shown in gray. 



2010) observations (A = 1135 - 1795 A) as well as 
archival data at 905 — 1187 A from the Far Ultraviolet 
Spectroscopic Explorer {FUSE; Moos et al. 2000; Sah- 
now et al. 2000). We confirm the featureless power- 
law nature of the spectrum over this wavelength range. 
Absorption is seen in 42 intervening systems including 
41 Lya absorbers and six metal-line systems. The fre- 
quency of IGM absorbers is consistent with larger sur- 
veys using FUSE&nA HST/STIS data (Danforth & Shuh 
2005, 2008; hereafter DS08), and the systems are spread 
across the entire redshift range covered by the combined 
COS/Ff/S'.B dataset (z < 0.47). 

The observations and data reduction techniques are 
discussed in §2, and we present a preliminary catalog of 
absorption lines in §3. Our conclusions are presented in 
§4. 

2. OBSERVATIONS AND DATA ANALYSIS 

Far-UV observations of lES 1553+113 were carried out 
2009 September 22 by HST/COS as part of the COS 
Guaranteed Time Observations (PID 11528, PI Green). 
Five exposures were made in each of the G130M (1135 < 
A < 1480 A) and G160M (1400 < A < 1795 A) medium- 
resolution gratings [R « 18,000) totalling 3.1 and 3.8 
ksec, respectively. Four central wavelength settings at 
each grating dithered known instrumental features along 
the spectrum and provided continuous spectral coverage 
over 1135 < A < 1795 A (see Green et al. 2010; Osterman 
et al. 2010). After retrieval from the archive, all ten 
exposures were reduced locally using CalCOS v2.11f. 

Flat-fielding, alignment, and coaddition of the pro- 
cessed exposures were carried out using IDL routines de- 
veloped by the COS GTO team specifically for COS FUV 
data^. First, the data were corrected for the most egre- 
gious instrumental features. While attempts at a true 
"flat-fielding" of COS data show promise, the technique 
is not yet robust enough to improve data of moderate 
S/N. However, we are able to correct the narrow ^15%- 
opaque features arising from ion repellor grid wires in 

^ Sec http://casa.colorado.edu/~danforth/costools.litml 
for our coaddition and fiat-fielding algorithm and additional dis- 



the detector. A one-dimensional map of grid-wire opac- 
ity for each detector was shifted from detector coordi- 
nates into wavelength space and divided from the flux 
and error vectors. Exposure time in the location of grid 
wires was decreased to ~ 70%, giving these pixels less 
weight in the final coaddition. We also modify the error 
and local exposure time at the edges of the detector seg- 
ments to de- weight fiux contributions from these regions. 
With four different central wavelength settings per grat- 
ing, any residual instrumental artifacts from grid-wire 
shadows and detector segment boundaries should have 
negligible effect on the flnal spectrum. 

The exposures are aligned with each other and interpo- 
lated onto a common wavelength scale. One exposure in 
each grating/detector was picked as a wavelength refer- 
ence, and the remaining exposures were cross-correlated 
with it. The wavelength region of cross-correlation for 
each case was picked to include a strong ISM absorption 
feature, and shifts were typically on the order of a reso- 
lution element (^ 0.07 A) or less. The COS wavelength 
solution has not yet been rigorously characterized, and 
we see a systematic shift between strong ISA4 lines and 
their expected LSR velocities. The shift is approximately 
constant across the COS wavelength range, so we apply 
a uniform +0.17 A shift to the wavelength vectors {^ 40 
km s^^ at ^ 1300 A) to bring ISM line centroids to the 
expected wlsr ~ seen in many ISM absorbers. 

Next, the aligned exposures were interpolated onto a 
uniform wavelength grid and coadded. The flux at each 
position was taken to be the exposure-weighted mean of 
flux in each exposure. Since exposure time was reduced 
in certain wavelength locations, as noted above, pixels 
near detector edges and where grid-wire shadows were 
removed received less weight than those in less suspect 
locations. The combined data show S/N ~ 20 per 7- 
pixel {^ 0.07 A) resolution element and are sufficient to 
detect narrow absorption features down to Wx « 15 mA 
at 4(7 significance. Figure 1 shows the entire combined 
COS/G130M and COS/G160M spectra. Figure 2 shows 
a more detailed view of the spectrum with prominent 
lines marked. 

In addition to the COS data, we utilize 45 ksec of Far 
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Fig. 2. — A more detailed view of the COS/G130M+G160M dataset smoothed by 7 pixels (approximately one resolution element). Error 
is shown in gray. Prominent ISM lines are marked with plus signs. IGM Lyo absorbers are marked with large vertical ticks. Smaller ticks 
denote corresponding Ly/3 detections. IGM metal absorbers are marked with open circles. The two question marks denote the ambiguous 
features discussed in §3. See Table 1 for line identifications and measurements. 
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Ultraviolet Spectroscopic Explorer (FUSE) observations 
taken 2004 April as part of program E526 (PI: Savage). 
While FUSE data alone are insufficient to characterize 
the H I absorber systems along a sight line, far-UV cov- 
erage is invaluable for confirming Lya lines at z < 0.11 
via Ly/3 absorption. Additionally, O VI AA1032, 1038 
and C III A977 absorbers are found only in FUSE data 
at z < 0.10 and z < 0.16, respectively. Thirty-seven 
FUSE exposures were retrieved from the archive and 
processed in the usual manner (Danforth & ShuU 2005). 
The final FUSE spectrum covers A = 905 - 1187 A with 
S/N « 7 — 10 per ^ 20 km s^^ resolution element. 

3. ABSORPTION LINES 

An initial analysis of the spectrum reveals a wealth 
of far-UV absorption features (Figures 1 and 2). Many 
of these are clearly Galactic ISM lines typical of most 
sight lines to Galactic and extragalactic sources. We 
label the remainder as redshifted IGM absorbers. To 
identify these lines, we follow a procedure similar to that 
employed in DS08: starting from the long- wavelength 
end of the spectrum, we interactively mark the strongest 
absorpton features, tentatively identifying them as Lya. 
The location of the corresponding Ly/? absorption is then 
checked, as are those of prominent metal-ion absorbers 
(O VI AA1032, 1038; C IV AA1548,1550; Si III A1207; 
C III A977, etc.). This process is iterative, as we iden- 
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tify weaker and weaker features. If there is component 
structure in a line profile that can be unambiguously de- 
convolved into multiple absorbers, we list these systems 
separately. However, most systems are listed as a single 
absorber, even if they possess rather complex line pro- 
files. 

We note two significant absorption features at 1197.25 
A and 1645.9 A with highly ambiguous identifications. 
The weak feature at 1197.25 A (T^^obs ~ 50 mA) cannot 
be Lya, nor is it consistent with either a higher-order 
Lyman line or any obvious metal-ion absorber for any of 
the known H I systems. The most plausible identifica- 
tion is that of O VI A1038 at z = 0.1538 (Fig. 3). The 
stronger A1032 line of the O VI doublet is blended with 
Galactic Si II A1190. No H I absorption is seen at this 
redshift in Ly/?, and a 4(t upper limit on the column- 
density can be set at log A^hi < 13.81. Lya absorption 
at this redshift is blended with the weaker line of the 
Galactic Si IV doublet at A1403. However, the Galactic 
Si IV lines appear in the expected 2:1 ratio, leaving little 
room for additional blended Lya z — 0.1538 absorption. 
It appears possible that this is a WHIM absorber with 
high enough temperature and metallicity that no neutral 
gas is seen (see also Savage et al. 2010). 

The strong absorption line at 1645.9 A (Wobs ~ 357 
mA) could be identified as Lya at z — 0.3539, but the 
expected Ly/3 absorber (Wobs > 50 mA) is not seen 
(Fig. 4a, b). The line is consistent with being Ly/3 ab- 
sorption a,t z — 0.6046, and an equivalent Ly7 feature is 
seen at 1560.3 A at approximately the expected strength 
(Fig. 4a, c). However, the latter feature is consistent with 
Galactic C I absorption lines seen elsewhere in the data 
(Fig. 4d-f). Therefore, we tentatively identify this fea- 
ture as a multi-component Lya system at z = 0.3539. 
Two or more Ly/? features of the required strengths can 
plausibly be hidden in the noise at the required location. 

Table 1 lists measurements for all detected IGM ab- 
sorption lines, grouped by redshift, and including the two 
ambiguous cases above. In total, we identify 42 IGM ab- 
sorbers (Table 1), 41 of which are detected in at least 
Lya. Corresponding higher-order Lyman line and/or 
metal ion absorption is seen in 15 absorbers. Seven sys- 
tems show metal absorption. The observed Lya absorber 
frequency per unit redshift, df\f/dz « 87 ± 15, down to 
a limiting equivalent width of 50 mA (~ 10"'^^ cm^^), is 
similar to that found for the larger DS08 sample to the 
same limit {dAf/dz = 95 ± 5). A more thorough search 
for broad Lya absorbers with 5 > 40 km s^^ will be con- 
ducted, following the receipt of additional data on this 
source scheduled for Cycle 18. Therefore, we caution the 
reader that this line list may not be complete for lines 
with 6 > 40 km s^^. 



OVI 1038 
z=0.153B 




logN„„=13.85_ 
bovi=26 km/s_ 


1196.5 


1197.0 1197.5 
Wavelength (A) 


1198.0 



0.0 



Fig. 3.— A weak absorption line at 1197.25 A {Wx = 50±9 mA) 
is consistent with OVI A1038 at 2 = 0.1538 (bottom panel). The 
corresponding O VI A1032 line is blended with Galactic Si II+S III 
absorption, but consistent with the measurement from the O VI 
A1038 line (middle panel). The corresponding Lya absorber is 
blended with Galactic Si IV A1403, but no H I is seen in Ly/3 (top 
panel) to 4a limit log A^hi < 13.81. 



4. RESULTS AND DISCUSSION 

An additional six orbits of COS integration time 
planned for Cycle 18 should improve the S/N of the com- 
bined dataset by a factor of ^ 2. Greatly improved S/N, 
as well as our evolving understanding of the COS instru- 
mental effects, will enable us to reliably measure low- 
contrast absorbers such as broad Lya systems and weak 
metal lines. We defer a more exhaustive analysis of the 
sight line until then, but note two key results here. 
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Fig. 4. — A broad feature at 1645.9 A (panel a) is interpreted as Lya at 2 = 0.3539. However, the expected Ly/3 feature (panel b) is 
not seen in the data. If the 1645.9 A feature is instead interpreted as Ly/3 at z = 0.6046, a feature of approximately the correct strength 
is seen at 1560.3 A, the location of the expected Ly7 at z = 0.6046 (panel c). However, the 1560.3 A feature is consistent with Galactic 
CI absorption (dashed profile in panels d-f). We conclude that the 1645.9 A feature is most likely a multi-component Lyo absorber at 
z ai 0.3539, but the individual components are too weak to appear in Ly/3. 



4.1. Triple Absorber Complex at z — 0.188 

The most interesting of the previously undiscovered 
IGM absorption systems is the triplet of metal-rich ab- 
sorbers at z = 0.18640, 0.18773, and 0.18989 (Figure 5). 
In Lya, the strong central absorber a.t z — 0.18773 is 
flanked by two weaker components at z — 0.18640 and 
z = 0.18989, orv = -399 km s'^ and v = +648 km s^^, 
respectively, relative to the system at z = 0.18773. 

The three absorption systems at z ^ 0.188 span ^ 1000 
km s~^ in comoving velocity space, appropriate for a 
large-scale filament in the cosmic web. We searched the 
SDSS Data Release 7 galaxy redshift catalog (Abaza- 
jian et al. 2009) for galaxies within one degree of the 
lES 1553+113 sight line and plotted them as a func- 
tion of redshift (Fig. 6, top panel). While the SDSS 
is complete only in the brightest galaxies (L > 3 L*) 
at this redshift, a clear concentration appears at z = 
0.187+ 0.003. None of these galaxies is closer than 24' 
{^ 4.5 Mpc at z = 0.188) to the fine of sight, so it 
is hard to claim a specific galaxy-absorber relationship 
(Figure 6, bottom panel). However, the median redshift 
of the galaxy sample is z = 0.187 and the la deviation 
(ctz = 0.0027) is roughly one-third of the redshift search 
space (Az = +0.008). This tight clustering around the 
absorber redshift, as well as the observed spatial distri- 
bution of the brightest galaxies, suggests that the galax- 
ies trace a large-scale filament in the cosmic web and 
that the absorption in the COS observations arises in 
the same structure. Deeper galaxy survey work (Keeney 
et al., in preparation) is complete to much lower luminos- 
ity {L > 0.3 L*) and may show a closer galaxy-absorber 
relationship. 



O VI absorption is seen in all three systems 
(log Novi = 13.4 + 0.3, 14.1 + 0.1, and 13.5 + 0.1, re- 
spectively). Strong N V absorption is seen in the cen- 
tral component (log TVnv = 13.7 + 0.1). DS08 measure 
(Nnv/Novi) = 0.24+°;22 in eleven O VI+N V low-z IGM 
absorbers, so the ratio observed at z ~ 0.18773 toward 
lES 1553+113 (Nnv/Novi == 0.4) is high but within the 
observed range. The DS08 sample is a biased toward 
higher TYnv values, the observed ratio suggests an ele- 
vated N/0 abundance in this absorber. 

C III is detected in the central and blue components 
(log Ncm = 13.25+2;o4 and 12.64 + 0.15, respectively), 
but 4(7 upper limits of log iVsan < 11.64 and log Nsny < 
12.29 can be placed on Si-ion absorption in all three sys- 
tems. Si II A1260 is tentatively detected at z = 0.1877 
as a pair of weak, narrow components with a total col- 
umn density log iVsin ~ 12.1. However, we do not de- 
tect other Si II lines, nor equivalent absorption in C II 
A1334.5 or C II A1036.3 (log Nqii < 12.82) and other 
singly-ionized species. 

It is likely that the gas in the central z ~ 0.1877 Lya 
system is multi-phase in nature, with a warm-hot ionized 
medium (WHIM) component traced by O VI and N V 
and a cooler, photoionized component traced by H I and 
C III. The "multiphase ratios" for these absorbers (Dan- 
forth & ShuU 2005, 2008) are Nm/Novi -- 1.3, - 0.6, 
and ^ 0.5 for the three main components. Typical val- 
ues for absorbers with similar Nm are ^ 0.6, ^ 2.5, and 
~ 0.8, respectively (Danforth & ShuU 2005). We can use 
the Lya and low-ionization metal detections and upper 
limits to constrain metallicity and relative abundances 
in the photoionized gas. In particular, C III and Si III 
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Fig. 5. — Detail of the absorption complex at 2: = 0.188. Top left panel shows three strong Lya absorption systems at z = 0.1864, 
z = 0.1877, and 2 = 0.1897. The reddest of these is clearly split into two components (2 = 0.18958 and 2 = 0.18989). The strong, 
central absorber also shows clear evidence of multiple structure, but is harder to deconvolve unambiguously. In subsequent panels, these 
components are marked by arrows. Corresponding metal absorption is detected in several of these components in C III, O VI, and N V, 
however there is no measurable absorption in either Si 111 or Si IV. Several unrelated lines appear in the profiles and arc labeled with vertical 
ticks. Data is binned by four pixels (~50% of a resolution element). 



have similar ionization potentials and are often detected 
in the same systems. At solar abundance ratios (As- 
plund et al. 2009), carbon is 8.3 times more abundant 
than silicon, but Si III is detectable to much lower col- 
umn densities than C III owing to the very strong /-value 
of the 1206.5 A transition (Shull et al. 2009). Thus, the 
two ions are often seen together in photoionized IGM 
absorbers (DS08). 

We measure iVciii/^Siiii > 40 in the z = 0.1877 ab- 
sorber, an unusually high lower limit. DS08 report C III 
and Si III detections in 22 low-z IGM systems with a 
median distribution of A^ciii/-^Siiii = 8.5lg°5^ (lo")- In 
Galactic High Velocity Clouds (Fox, Savage, & Wakker 
2006; ShuU et al. 2009) the ratio, A^ciii/A^Siiii, typically 
ranges from 5-20. Thus, the abnormally high ratio ^ 40 
found in these IGM absorbers is well outside the usual 
range. Comparing our measurements with a grid of sim- 
ple CLOUDY photoionization models (detailed in Dan- 
forth et al. 2006), we see that the relative column densi- 
ties of C III and Si III are fairly insensitive to photoion- 
ization parameter U = n^/nu- Typical model ratios 
are (A^ciii/-^Siiii) ~ 10 in the expected range of IGM 
photoionization parameters {U ~ 10~^) and are largely 
insensitive to assumptions about metallicity, photon con- 
tinuum, and gas density. 

The unusually high C Ill/Si III ratio suggests that the 
C/Si abundance in this system may have a strongly non- 
solar abundance pattern. If A*"siiii is typically ^10% of 



A^ciii: as seen in other IGM observations and models, 
[Si/C] > 0.6, or greater than four times the solar value. 
Comparing the observed Nciii/Nyii ratio with the mod- 
els, we expect log (A^^ciii/^hi) < -1-7 for Z = 0.1 Z©. 
However, the observed column density ratio is an order of 
magnitude higher, suggesting that (C/H) is close to solar 
values in this system. Without additional low-ionization 
species detected, we cannot determine whether carbon is 
overabundant or silicon is underabundant relative to the 
solar ratio. The C III detection at 2; = 0.1864 is factor of 
four weaker than in the main absorber, while the upper 
limit on A^siiii is the same. Therefore, this system puts 
weak constraints on the metallicity and abundances in 
this absorber. 

Ahhough (C/Si) ^ 8.3 ± 1.2 in the Sun (Asplund 
et al. 2009), variations in this abundance ratio can oc- 
cur, depending on the youth of the stellar population 
and its initial mass function (IMF). Carbon is produced 
primarily by helium burning in intermediate-mass stars 
(red giants, horizontal branch), whereas silicon arises 
from more advanced a-process nucleosynthesis in mas- 
sive stars. The usual abundance trends show enhanced 
(Si/C) and reduced (N/0) in low-metallicity stellar pop- 
ulations (McWiUiam 1997; Cayrel et al. 2004). Theo- 
retical predictions (Woosley & Weaver 1995) show that 
[a/Fe] increases with increasing progenitor mass (here, 
a includes O, Mg, Si, S, Ca, Ti). Thus, a low Si and 
O abundance compared to C and N suggests an IMF 
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Fig. 6. — Correlating the triple absorption system at z ~ 0.188 
with SDSS galaxies within one-degree of lES 1553+113. The top 
panel shows a the redshift distribution of SDSS galaxies within 
the search radius. The three Lyo absorber redshifts (arrows) cor- 
respond to a significant peak in the galaxy distribution bounded 
by the vertical dotted Unes (0.180 < z < 0.196). We plot these 
eight galaxies (diamonds) in relation to the lES 1553-1-113 sight 
line (star) in the lower panel. Spectroscopic redshifts are noted 
next to each diamond, and symbol size denotes galaxy luminosity. 
The distribution of bright galaxies at approximately the redshift 
of the triple absorption system (z ^ 0.188) suggests a large-scale 
filament. Though none of the SDSS galaxies are closer than 24' 
(^ 4.5 Mpc at z = 0.188) to the Une of sight, the SDSS galaxy cat- 
alog is not complete to lower-luminosity objects at this redshift. 
All galaxies in this field have L > 3 L* . A deep galaxy redshift 
survey will likely show fainter galaxies closer to the AGN line of 
sight. 



skewed toward low-mass stars. 

Comparing H I or C III to high ions O VI and N V 
requires them to be in the same thermal phase for a 
meaningful analysis. Hybrid ionization modeling (CIE 
plus photoionization) of the high ions in this system are 
reported by Yao et al. (2010), in which the z = 0.1877 
system is used as a test case for a physical parameter- 
based absorption line modeling exercise. 

4.2. Constraining the Redshift of lES 1553+113 

The redshift of lES 1553-1-113 is crucial to determining 
the intrinsic properties of the source. Indirect methods 
of constraining the redshift of lES 1553+113 fall into two 
categories. First, the ratio of AGN to host galaxy optical 
luminosity in BLLacs is thought to cover a fairly small 
range (Wurtz, Stocke & Yee 1996; Sbarufatti, Treves, 
& Falomo 2005). Various deep ground-based (Hutch- 
ings & Neff 1992; Scarpa et al. 2000) and space-based 
(Urry et al. 2000; Carangelo et al. 2003; Sbarufatti et al. 
2006; Treves, Falomo, & Uslenghi 2007) optical studies 
have failed to detect a host galaxy beneath the glare 
of the AGN. From these non-detections, redshift lim- 



its from z > 0.09 to z > 0.78 have been set by vari- 
ous groups. Treves, Falomo, & Uslenghi (2007) refined 
this tO2; = 0.3 — 0.4 using a more sophisticated analysis 
of the same optical data. However, the validity of the 
assumption of host/nuclear luminosity relationship has 
been called into question by O'Dowd & Urry (2005). 

A complementary technique uses the observed very 
high energy (VHE) spectrum (0.1 — 10'^ GeV) to place 
limits on the redshift of BL Lacs. This method assumes 
that the VHE spectral energy distribution (SED) of an 
object will be modified, as TeV photons interact with 
photons in the ambient extragalactic background and 
produce e+e^ pairs (e.g.. Younger & Hopkins 2010; Per- 
sic & De Angelis 2008). The longer the pathlength, the 
steeper the VHE SED becomes. Uncertainties in the 
extragalactic IR background and the intrinsic SED of 
the AGN render this method uncertain, but the red- 
shift of lES 1553+113 has variously been constrained to 
z < 0.74 (Aharonian et al. 2006; Albert et al. 2007) or 
z < 0.80 or z < 0.42 (Mazin & Goebel 2007) based 
on HESS and MAGIC observations. Abdo et al. (2010) 
use data from the first six months of Fermi 7-ray ob- 
servations in conjunction with observations from radio 
wavelengths to 1 TeV to model the intrinsic SED of 
lES 1553+113. Based on these models, they determine 
a redshift z — 0.751qq5. The error bars on this estimate 
appear to be much smaller than justified by this method. 

From Figures 1 and 2, it is clear that there are H I 
systems throughout the redshift range from z ^ to 
near the end of the COS spectral coverage {z = 0.47). A 
strong line at 1695 A is identified as Lya at z = 0.395 
and confirmed by detection of O VI at the same redshift 
in both lines of the doublet. This sets a firm lower limit 
on Zcm > 0.395. Two weaker features at 1709.5 A and 
1741.5 A appear in the data, which we identify as Lya at 
z = 0.4063 and z — 0.4326, respectively. Though we do 
not detect higher-order Lyman or metal-ion lines at these 
redshifts, the two z > 0.4 Lya absorbers are weak enough 
that we do not expect confirmation in other lines. The 
continuum of the BL Lac object remains smooth across 
the entire COS band (Figs. 1 and 2), and no intrinsic 
emission or absorption is seen. 

Thus, we can confidently constrain the emission red- 
shift of lES 1553+113 to Zom > 0.400, and it appears 
likely that it may be as high as Zem = 0.433. Addi- 
tional COS observations may detect corresponding Ly/? 
absorption to these two z > 0.4 absorbers (WLy/3 ~ 12 
mA is expected). The confirmed and unconfirmed di- 
rect redshift limits from the HST/COS observations are 
compatible with both the lower limits set by the non- 
detection of an optical host galaxy (zcm ^ 0.1 — 0.4; 
Urry et al. 2000; Sbarufatti et al. 2006; Treves, Falomo, 
& Uslenghi 2007) and the VHE SED upper limits upper 
limits {zem ^ 0.8; Aharonian et al. 2006; Albert et al. 
2007; Mazin & Goebel 2007). 

We now assess the validity of the most recent VHE 
SED redshift estimate (zcm = O.75+005) from Abdo ct al. 
(2010). Our current COS far-UV spectra (G130M and 
G160M) are only sensitive to Lya absorbers at z < 0.47. 
However, higher redshift absorbers can be detected us- 
ing the less sensitive higher-order Lyman lines or O VI 
doublet. COS far-UV data cover the wavelength range 
1135 A < A <1795 A corresponding to Ly/3 redshifts 
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Fig. 7. — Observed and expected distribution of IGM absorbers along the lES 1553+113 sight line. Vertical ticks mark the redshift of 
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be detected in the COS far-UV band and we must rely on much less sensitive Ly/3 detections. As discussed in the text, no z > 0.4 Ly/3 
absorbers are detected. 



0.11 <z< 0.75 and O VI redshifts 0.10 <z< 0.74. Ly/3 
and O VI systems at z > 0.47 will appear at A > 1508 A. 
We find empirically (DS08) that the detection threshhold 
for absorption lines with no prior "signposts" (such as 
known absorber redshifts) is - lOcr. The COS/G160M 
data in this region are of sufficiently high quality that 
we would expect to detect lines of PFobs ^ 50 mA at 
a ^ lOcr significance level. This corresponds to rest- 
frame Wr > 30 mA for Ly/3 absorbers at z ^ 0.5 or 
log ^Hi > 13.6 for Ly/3. 

Figure 7 shows the observed and predicted distribution 
of IGM absorbers as a function of redshift. We calculate 
the expected number of absorbers A/'abs per Az = 0.025 
bin (alternately, dAf/dz) based on the S/N-determined 
minimum equivalent width in each bin and the sample of 
^ 650 H I absorbers from DS08. The dashed curve in the 
figure represents no dM/dz evolution with redshift. The 
evolution of low-z H I absorbers is somewhat uncertain. 
However, if we assume that the H I absorber frequency 
evolves as {dJ\f/dz)iii oc (l + z)'' (Penton, Stocke, & ShuU 
2004) and adopt 7^0.7 for a modest evolution between 
z = and z ^ 1, the expected number of H I detec- 
tions (dotted curve in Figure 7) rises at higher redshift 
by ^ 20 — 50%. The sharp drop in expected detections 
at z > 0.47 coincides with the switch from Lya to Ly/3 
as an IGM tracer (shaded area) and the resulting loss 
of sensitivity discussed above. Summing the expected 
number of Ly/3 absorbers (A/'abs.exp) which should ap- 
pear at A > 1508 A data from Figure 7 over the range 
0.47 < Zabs < 0.75, we find 7Vabs,cxp ~ 7 and Wabs.oxp ~ 9 
in the constant and evolved H I models, respectively. 



Thus, we should expect ^ 8 high-z Ly/3 absorbers in the 
lES 1553+113 data if Zom > 0.75. 

Are the predicted high-z Ly/3 absorbers seen? As dis- 
cussed in §3, the strong feature at 1645.9 A can be ruled 
out as Ly/3 at z — 0.6046 since the corresponding Ly7 line 
is consistent with Galactic C I (Fig. 4). Thirteen other 
absorption features longward of 1507 A have been identi- 
fied as Lya lines, and eight of these identifications are not 
confirmed with higher-order Lyman or metal-ion lines. If 
these eight single-line detections are instead interpreted 
as potential high-z Ly/3 systems, six have inconsistent 
Ly7 non-detections and a seventh has Ly7 blended with 
another line. The only possible high-z Ly/3 absorber is a 
marginal line detected at 1584.3 A and identified as Lya 
at z = 0.30328. If this line is instead identified as Ly/3 
at z — 0.54463, the corresponding Ly7 absorber should 
be at 1502.2 A, nearly concident with a weak feature 
identified as Lya at z = 0.23559. The relative strengths 
of the two features are consistent with H I absorbers at 
z = 0.54463, but both line detections are of relatively 
low significance. Additional COS observations may im- 
prove the significance of the line detections. However, 
we find none of the other predicted high-z H I or O VI 
absorbers in the data, so the Zem > 0.545 redshift limit 
from the possible Ly/3 detection for lES 1553+113 is very 
speculative. 

The additional far-UV observations of lES 1553+113 
scheduled for HST Cycle 18 will improve the S/N of 
the current dataset by a factor of ~ 1.7 and conse- 
quently lower the minimum detectable line strength for 
potential Ly/3 absorbers at z > 0.47 by a similar fac- 
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tor. However, observations with the COS/G185M grat- 
ing covering the wavelength range 1800 A < A < 2100 A 
(0.47 < Zabs S 0-'''3 in Lya) would be more efficient 
at detecting z > 0.47 IGM absorbers in the sight line. 
Despite the relatively lower efficiency of the COS near- 
UV gratings and detectors compared with their far-UV 
counterparts, the Lya lines will be ^ 3 — 7 times stronger 
than Ly/3 and should be easily detected or ruled out with 
only a few kiloseconds of observations. Such observations 
could potentially also measure weak intrinsic broad-line 
Lya emission from the BL Lac object, as has been seen in 
other BL Lac objects at lower redshift (Danforth, Stocke, 
Winter, et al., in prep). 

We constrain the source redshift of lES 1553+113 sta- 
tistically by truncating the A/'abs, expected rnodel curves 
in Figure 7 at a range of 0.4 < Zcm < 0.75. Applying 
a Kolmogorov-Sniirnov test to the different models, we 
set a Icr constraint of Zcm < 0.58 for the non-evolved 
model (and Zom S 0-49 for the evolved H I distribution). 
An emission redshift of ^ 0.75 is ruled out for both 
models at a 90% or greater level of confidence. Thus, 
we constrain the redshift of lES 1553+113 to the range 
0.43 < Zem < 0.58. 

5. CONCLUSIONS 

The BLLac object lES 1553+113 is one of the bright- 
est objects in the sky in 7-rays, as well as being a no- 
table UV and X-ray source. However, the AGN emis- 
sion is that of a relativistic jet aligned closely with our 
line of sight and, like most such objects, has no intrinsic 
emission or absorption features at any wavelength. This 
featureless, power-law continuum is ideal for measuring 
intervening IGM features that are weak and broad, such 
as thermally-broadened Lya systems. However, the lack 
of intrinsic features makes constraining the redshift of 
the object difficult. 

We present unprecedented high-quality far-UV 
HST/COS and FUSE spectra of the BLLac object 
lES 1553+113 at spectral resolution 15-20 km s~^. 
These data show 42 intervening IGM absorbers, 41 of 
which are detected in Lya, and 15 in Ly/3 and/or metal 
lines. The richest absorption system in the line of sight 
is a trio of Lya absorbers at 2; « 0.188 covering ^ 1000 
km s~^ of velocity space. Several metal ions are also 
detected in these systems, including O VI, N V, and 
C III. However, neither Si IV nor Si III is detected in any 
of the systems. The C Ill/Si III ratio implies a (C/Si) 
abundance at least four times the solar value, while a 
high N V/0 VI value suggests an overabundance of N 
as well. A detailed analysis of the physical conditions in 
this system can be found in Yao et al. (2010). 

The redshift of lES 1553+113 has never been deter- 
mined directly, and the only limits placed on it come 
from indirect means such as the shape of the 7-ray spec- 
trum and the lack of an AGN host galaxy in deep optical 
images. A strong Lya+0 VI absorber at 2; = 0.3951 
gives the first direct lower limit to the redshift of the 
object. Two weaker Lya absorbers at z = 0.4063 and 
z = 0.4326 give slightly higher estimates of the redshift, 
but these weak Lya lines are not confirmed by additional 
line detections. 

These lower limits are consistent with most previous 
measurements via optical non-detections of host galax- 



ies and 7-ray SED constraints. Abdo et al. (2010) de- 
rive Zcm — 0.75to;o5 based on the latest Fermi and TeV 
7-ray SED measurements, considerably higher than our 
intervening absorber upper limits. COS far-UV spectra 
are not sensitive to Lya absorbers at z > 0.47, but the 
G160M grating has some sensitivity to intervening Ly/3 
and O VI absorbers out to z ^ 0.75. If the Abdo et al. 
(2010) redshift estimate were accurate, we would expect 
to find -- 8 Ly/3 absorbers at 0.47 < z < 0.75. We find 
no evidence for any higher redshift absorption systems. 
There are only a few absorption features at A > 1500 A 
with ambiguous line identifications that could potentially 
be Ly/3 systems at z > 0.47. While these systems are in- 
dividually suggestive, we find nowhere near the number 
of absorbers predicted statistically. We conclude that 
the redshift of lES 1553+113 is not much higher than 
z « 0.45. 

lES 1553+113 is one of the brightest X-ray sources on 
the sky and has been suggested as a sight line that could 
be efficiently probed for WHIM absorption in O VII. The 
combined O VI column density in the three absorbers at 
z ^ 0.19 is ^ 2 X lO^** cm~^. Spectrographs on modern 
X-ray observatories are sensitive to log A^ovii ^ 15.5. If 
the temperature of any of these O VI systems is high 
enough, sufficiently long Chandra and/or XMM/Newton 
observations may reveal a O VII counterpart to these 
O VI absorbers that could constrain the long-sought X- 
ray WHIM (Bregman 2007). However, at the observed 
Li-like (O VI) oxygen column density, log A"ovi ~ 14.3 
in the trio of absorbers, the expected column densities of 
He-like (O VII) and H-like (O VIII) oxygen are probably 
just below the detectability levels of Chandra and XMM. 
Recent analysis of stacked X-ray absorption data (Yao 
et al. 2009) at the known IGM redshifts of O VI absorbers 
finds no evidence for O VII or O VIII absorbers to a limit 
A^oviiZ-^ovi < 10. Therefore, the z « 0.19 absorbers 
might have X-ray column densities log iVovii < 15.3, 
just below the limits of current X-ray observatories. 

Finally, these observations showcase the powerful new 
tool available to astronomers for probing the low-redshift 
IGM. COS is 10 — 20 times more sensitive in the far-UV 
to point sources than previous instruments on HST. An 
additional six orbits of COS observations arc planned 
for 2010, which should improve the S/N of the combined 
dataset by a factor of ^ yS- Improving the data qual- 
ity will help confirm or refute some of the tentative line 
identifications from this paper and will undoubtedly un- 
cover additional weak absorbers. We will place further 
constraints on [C/Si] and [N/0] in the z = 0.188 system, 
identify new broad, shallow Lya absorbers, and investi- 
gate possible high-z Ly/3 systems with our new Cycle 18 
observations. 
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of people who made the HST Servicing Mission 4 the 
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Penton, Stephane Beland, and the other members of the 
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calibration and verification. C. D. wishes to acknowl- 
edge a fruitful discussion with members of the KIPAC 
consortium. This work was supported by NASA grants 
NNX08AC146 and NAS5-98043 to the University of Col- 
orado at Boulder. 
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TABLE 1 

Intergalactic Line Detections 



2abs 


Line 


Wr (mA) 


b (kms-l) 


log Af 


notes 


0.00717 


Lya 


108 ± 12 


35 ±4 


13.37 ±0.04 


blend with Viz = 0.187 


0.02386 


Lya 


54 ±8 


19 ±4 


13.07 ±0.05 




0.03254 


Lya 


24 ±8 


13 ±7 


12.70 ±0.10 


weak 


0.04012 


Lya 


215 ±9 


27 ±1 


13.86 ±0.02 






Ly/S'' 


87 ±46 


11±6 


14.39 ±0.22 


ambiguous ID 








17+f 


14.25«-I? 


H I COG solution'' 


0.04281 


Lya 


135 ± 14 


73 ±'7 


13.44 ± 0.04 


broad 


0.04498 


Lya 


275 ± 10 


47 ±2 


13.87 ±0.01 


asymmetric 




Ly/3- 


90 ±8 


45 ±2 


14 iq+0.09 










24tl 


14.20 ±0.07 


H I COG solution'' 


0.05094 


Lya 


84 ±9 


17 ±3 


13.32 ±0.04 




0.07832 


Lya 


58 ±10 


24 ±6 


13.09 ±0.06 


weak 


0.09481 


Lya 


176 ±9 


30 ±2 


13.67 ±0.02 






Ly/3- 


34 ±8 


26 ±11 


13.8010-1^ 


weak, ambiguous ID 








2StTo 


13.70+111 


HI COG solution'' 


0.10230 


Lya 


360 ±4 


52 ±1 


14.12lH^ 


BLA? 




Ly/3- 


81 ±10 


36 ±12 


14.2010:1? 










421- 


14.10 ±0.08 


H I COG solution'' 


0.12040 


Lya 


32 ±7 


21 ±9 


12.85ini 


weak 


0.12325 


Lya 


77 ±17 


81 ±16 


13.18 ±0.07 


weak, broad 




Si III 1206 


16 ±8 


9±8 


11.93 ±0.14 


tentative 


0.13334 


Lya 


70 ±13 


46 ±8 


13.15 ±0.06 




0.13483 


Lya 


156 ±9 


30 ±2 


13.60 ±0.02 




0.14028 


Lya 


28 ±8 


18 ±8 


12.76 ±0.10 


weak 


0.15164 


Lya 


51 ±9 


28 ±6 


13.02 ±0.06 




0.15234 


Lya 


173 ± 13 


66 ±5 


13.56 ±0.03 




0.15380 


O VI 1038 


50 ±9 


26 ±6 


13.95 ±0.06 


tentative, see §3, Fig. 3 


0.16540 


Lya 


35 ±8 


40 ±11 


12.83 ±0.09 


weak 
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Intergalactic Line Detections 



0.18640 


Lya 


139 ±7 


32 ±2 


13.52 ±0.02 






Ly7 


10 ±7 


~5 


13.84 ±0.56 


weak 








17+- 


i3.67+;;-^7 


poorly-constrained HI COG*^ 




O VI 1038 


10 ±4 


~5 


13.40 ±0.24 


low significance 




C III 977 


26 ±11 


22 ±14 


12.64 ±0.15 


noisy 


0.18773 


Lya 


293 ±8 


50 ±1 


13.89 ±0.01 


double 




Ly7 


23 ±32 


15 ±12 


14.02 ±0.27 


marginal 








361^4 


13.99lO:f, 


poorly-constrained HI COG'' 




O VI 1032 


108 ± 16 


41 ±6 


14.04«;»« 


complicated structure 




O VI 1038 


68 ±2 


35 ±1 


14.11«:«1 






NV 1238 


74 ±7 


36 ±4 


13.62«:«I 


two components 




N V 1242 


50 ±11 


36 ±9 


13.71 ±0.08 


single component 




cm 977 


89 ±6 


36 ±5 


10 oc+o.oe 

13-^O_0.04 


two components 




Sill 1260 


~ 18 


~5 


12.1 ±0.2 


two components, tentative 


0.18958 


Lya 


85 ±17 


49 ± 13 


13.24 ±0.11 




0.18989 


Lya 


74 ±8 


28 ±4 


13.20 ±0.05 






Ly7 


21 ±10 


7 ±10 


14.03 ±0.15 


weak 




OVI 1032 


38 ±8 


21 ±6 


13.53 ±0.07 






O VI 1038 


16 ±5 


13 ±9 


13.48tri9 




0.19503 


Lya 


124 ±7 


32 ±2 


13.45 ±0.02 




0.20621 


Lya 


58 ±12 


44 ±10 


13.06 ±0.07 


weak 


0.20747 


Lya 


53 ±12 


44 ±11 


13.02 ±0.08 


weak 


0.21640 


Lya 


87 ±14 


46 ±7 


13.25 ±0.05 




0.23559 


Lya 


20 ±7 


12 ±8 


12.62 ±0.12 


weak 


0.24155 


Lya 


242 ± 10 


36 ±2 


13.84 ±0.02 


strong, multi-component? 




Ly/3 


34 ±5 


26 ±7 


10 70+0.07 


weak 


0.25439 


Lya 


68 ±2 


31 ±2 


10 17+0.04 

^■^■^1 -0.05 




0.26560 


Lya 


46 ±10 


26 ±7 


12.97 ±0.08 


weak 


0.27753 


Lya 


77 ±10 


26 ±4 


13.22 ±0.05 




0.28944 


Lya 


56 ±10 


24 ±6 


13.07 ±0.07 




0.29830 


Lya 


56 ±12 


19 ±5 


13.09 ±0.07 




0.30328 


Lya 


70 ±7 


34 ±6 


,0 nn+0.06 


weak, noisy 


0.32235 


Lya 


190 ±4 


34 ±4 


13.73 ±0.04 






Ly/3 


55 ±9 


38 ±7 


13.92 ±0.06 










I8I2 


13.98t»:i? 


H I COG solution'' 


0.35390 


Lya 


235 ± 21 


~79 


13.73 ±0.03 


tentative, see §3, Fig. 4 


0.35650 


Lya 


170 ± 14 


43 ±4 


13.60 ±0.03 






Ly/3 


31 ±9 


28 ±9 


13.66 ±0.10 










(50+00 
28_i2 


13.66tJ!:}t 


poorly-constrained HI COG'' 


0.37884 


Lya 


136 ± 19 


47 ±6 


13.47 ±0.06 






Ly/3 


23 ±5 


15 ±5 


13.55 ±0.08 










27+5^ 


13.54 ±0.12 


poorly-constrained HI COG'' 


0.39505 


Lya 


278 ± 18 


57 ±4 


13.84 ±0.02 


multiple components 




Ly/3 


38 ±1 


~54 


10 7f5+0.07 


broad, multi-component 




O VI 1032 


57 ±6 


31 ±4 


13.71 ±0.04 






O VI 1038 


27 ±4 


23 ±5 


13.69«;«| 




0.40630 


Lya 


59 ±5 


29 ±2 


,0 10+0.11 

iJ.iJ_g Qg 


weak 


0.43261 


Lya 


89 ±6 


32 ±2 


10 01 +0.07 

ld.di_„gg 


noisy 



Measurement based on FUSE data. 

Curve of growth 6hi, -^hi solution based on multiple Lyman lines. 



